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Abstract 

Long-lived colored particles with masses m ^ 200 GeV are allowed by cur- 
rent accelerator searches, and are predicted by a number of scenarios for 
physics beyond the standard model. We argue that such "heavy partons" 
effectively have a geometrical cross section (of order 10 mb) for annihilation 
at temperatures below the QCD deconfinement transition. The annihila- 
tion process involves the formation of an intermediate bound state of two 
heavy partons with large orbital angular momentum. The bound state sub- 
sequently decays by losing energy and angular momentum to photon or 
pion emission, followed by annihilation of the heavy partons. This decay 
occurs before nucleosynthesis for m <^ 10^^ GeV for electrically charged 
partons and m ^ TeV for electrically neutral partons. This implies that 
heavy parton lifetimes as long as 10^^ sec are allowed even for heavy partons 
with m ~ TeV decaying to photons or hadrons with significant branching 
fraction. 



Fig. 1. A schematic representation of a hadron made of a heavy parton 
and surrounded by QCD "brown muck." 

The phenomenology and cosmology of heavy (m ^ GeV) long-lived colored particles 
has received renewed attention recently because of the proposal of "split supersym- 
metry" [1]. Another possible motivation is having a long lived gluino [2] or squark 
[3] as the next-to lightest superpartner in weak scale supersymmetry. More gener- 
ally, it is a phenomenologicall interesting possibility that such particles could exist. 
For m > 200 GeV such particles are consistent with current collider bounds, and if 
m < 2 TeV they will be accesible at LHC [4j. 

This note concerns the cosmology of such particles, which we refer to generically 
as "heavy partons," since they are constitutents of exotic long-lived hadrons at low 
energies. This has been studied by a number of authors, and there is some controversy 
about the the correct relic abundance, even at the level of the order of magnitude 
[5]. The disagreement is over the extent to which the strong interactions enhance 
the annihilation cross section in the early universe over the perturbative value. At 
temperatures below the deconfinement temperature Tc ^ 180 MeV, the heavy partons 
are confined inside hadrons. Just as in heavy quark effective theory pj it is useful to 
picture these hadrons as consisting of a heavy parton surrounded by QCD "brown 
muck" with a radius of order i?had ~ GeV~^ (see Fig. 1). If the heavy parton is a 
color triplet, the brown muck will involve at least one light quark, while if the heavy 
parton is a color octet, it may involve just gluons. Our arguments will not depend 
on the details of the brown muck. 

We therefore consider two heavy hadrons in the early universe at temperatures 
below the QCD phase transition. It is clear that the strong interactions give a geo- 
metrical cross section (of order -R^ad) ^r the heavy hadrons to interact, but this only 
means that the brown muck of one hadron interacts with that of the other. In order 
for the partons to annihilate, the wavefunctions of the heavy partons themselves must 
overlap, so the direct annihilation cross section is proportional to m~^. 
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Fig. 2. The formation of a highly excited bound state. 

We will argue that the strong interactions nonetheless give rise to an effective 
geometrical cross section for annihilation at temperatures below the QCD phase tran- 
sition. The first stage of this process is illustrated in Fig. 2. This is a capture process 
in which two heavy hadrons interact to form a bound state, with the energy carried 
off e.g. by a pion. As we will see, this bound state has large orbital angular momen- 
tum {L ~ 10 for m ~ TeV) and therefore does not decay promptly. In the second 
stage, the bound state loses energy and angular momentum, e.g. by emitting pions 
and/or photons. Bound states with binding energy larger than the temperature T 
survive collisions with particles in the thermal bath, and eventually decay through 
annihilation. The net result is that the heavy partons effectively have an enhanced 
cross section for annihilation given by the capture cross section. (This is similar to 
the scenario originally described in Ref. [1] for stable gluinos.) The phenomonological 
implications of this result will be discussed at the end of this note. 

We now explain why the capture process is unsuppressed. The first point is the 
existence of the bound state. The potential between two heavy partons can be written 
schematically as a sum of a short-distance Coulomb interaction and an attractive 
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linear term representing the effects of confinement: 

v[r) ~ or. (1) 

r 

Here a ~ string tension, and C is a group theory factor that depends 

on the color representation of the heavy partons. The color Coulomb force is always 
attractive in the color singlet channel {i.e. C < 0). We also expect the long-range part 
of the potential to be attractive in the color singlet channel, since it is responsible 
for color confinement. The energy spectrum of the system of two heavy partons 
therefore looks as shown schematically in Fig. 3. The low-lying states are Coulombic, 
with energy splittings of order aqc^m ^ Ahad ~ GeV, while the states near the 
continuum threshold are dominated by the linear term. Spin-dependent interactions 
are suppressed by 1/m, so spin excitations are small and we do not expect them to 
play an important role in the process we are considering. The rotational excitations 
are very important for determining the properties of the intermediate bound state. 
The minimum radius for a given angular momentum is determined by the circular 
orbits. In the regime where the linear term dominates, we have 

2 \ 1/3 

(2) 




The largest angular momentum occurs for r ~ -Rhad- At larger radii, the string 
confining the heavy partons will break and the spectrum becomes a continuum of 
two-particle states]^ The largest angular momentum of a bound state is therefore 

This estimate is consistent with the fact that there are expected to be L = 3 states in 
the the T system below the B-B threshold [7J. This gives Lmax ~ 3(m/m;,)^/^ ~ 40 
for m ~ TeV. 

We now consider the process depicted in Fig. 2, taking place at temperatures 
T < Tc. We are interested in impact parameters < Aj^^d; so that the brown muck 
clouds overlap, and therefore interact strongly. (We neglect the possibility that pion 
exchange gives a longer range attractive interaction.) To understand the possible 
transition to a bound state, we must understand what quantities are conserved in 
this transition, taking into account the special kinematics of the situation. Unlike a 
weakly interacting relic, a colored particle will be in kinetic equilibrium with gluons 

^String breaking is suppressed for large TVc- We will not keep track of factors of iVc in this paper. 
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Fig. 3. Spectrum of the system consisting of two heavy partons. (Spin 
excitations are not shown.) Bound states with large angular momentum 
will decay through a series of AL = 1 transitions, as shown. 



at the deconfinement temperature. The initial velocity is therefore determined by the 
temperature to be 



B 

m 



(4) 



where Tg is the temperature where the bound state is formed. Note that the momen- 
tum is large compared to the inverse range of the interaction 



Pi ~ mvi ~ [mTBY^'^ ^ 
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(5) 



had 



SO the scattering process is not dominated by s-wave scattering, and we expect higher 
partial waves to be important. On the other hand, the velocity is slow enough that 
the forces exerted by the brown muck can change the velocity significantly: 



At; aAt {XlJm){R^^Alvi) A 



had 



> 1. 



B 



(6) 
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We see that the velocity is not conserved in the coUision, even though the partons 
are heavy. 

Energy and angular momentum must of course be conserved in the collision. En- 
ergy conservation is satisfied by the emission of a pion (or perhaps several) that carries 
away the binding energy. The pion can also carry away some angular momentum, 
but only AL ~ few. The typical initial angular momentum is 

/ ™ X 1/2 /Try \ 

L,~m^,r,~(mT^)V2i?,^d~ 10 (^) (^j . (7) 

The initial angular momentum is large, but note that becomes smaller than L^^x 
for Tb <Tc. This means that the binding energy is 

/^27-2 \ 1/3 / 2r2\l/3 /^27-2 \ V3 



B=£.„-B,~^^j -y-^) ~A... (8) 

The initial kinetic energy is only ~ m^, but the binding energy is sufficiently large 
to produce a pion. 

It is also easy to see that the transverse distance between the heavy partons is 
allowed to change significantly during the collision. The forces of the brown muck 
give 

Ar^ ~ a^At^ ^ ^ f^)' i?had, (9) 

m \ V J Tb 

where ± refers to the component perpendicular to the line connecting the heavy 
partons. Conservation of angular momentum gives 

^ > 1. 10 

r V 

Based on the considerations above, we argue that the cross section for the forma- 
tion of the bound state is geometrical, i. e. 

The reason is simply that there is no symmetry or kinematic factor that suppresses 
the transition, so it should proceed at the going rate for strong interactions. In 
fact, it is possible that longer range forces due to pion exchange increase the cross 
section, a possibility that we will neglect here. As we have seen, the transition can 
occur for impact parameters of order i?had without violating any exact or approximate 
conservation law. Furthermore, the energy of the emitted pion(s) is of order Ahad; so 
there is no low-energy suppression. 
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It is important to keep in mind that the geometrical cross section gives the rate for 
the bound state to form, and we must consider the subsequent evolution of the bound 
state to see whether it enhances the annihilation of the heavy partons themselves. 
The first issue is whether bound states of the kind we are discussing can survive in the 
thermal bath long enough so that it can radiate away energy and angular momentum 
and eventually annihilate. The biggest threat to their survival is collisions with 
photons, since pions have low number density for T < ~ Tc. States with binding 
energy B > T cannot be destroyed efficiently, since the probability for finding a 
photon in the thermal bath with energy B is suppressed by e~^/^, which drops rapidly 
for T < B. (This rate will have additional suppression if the heavy partons and the 
brown muck are electrically neutral.) To find the precise value of the temperature 
below which destruction becomes inefficient requires a detailed hadronic model, but 
it is of order Tc. 

One might also worry about successive collisions with photons gradually increasing 
the energy of the bound state until it becomes unbound. The dominant process is 
inverse decay B^ B', where B and B' are bound states. However, these merely 
establish an equlibrium distribution of excitations (in which bound states are more 
numerous) on a time scale given by the rate for the decay B' Bj. 

We now consider the decay of the bound state. We consider first the case of elec- 
trically charged heavy partons, which can decay by the familiar process of photon 
radiation. The bound state decays by first radiating away energy and angular mo- 
mentum to get to a bound state with L ~ 1, which then decays by annihilation into 
quarks or gluons. We first estimate the decay rate to an L ~ 1 state. The hnear term 
dominates the force on a heavy parton for r ^ Rc, where 



. (12) 



This corresponds to a binding energy 



~ ^ ~ («QCD<7)^/^. (13) 

The time to lose energy of order Ec can be estimated by using the Larmor formula 
E ~ aa^, where a ~ a/m is the acceleration in the linear potential. Since the 
acceleration is constant, we have 

t(A£; ~ -E^ ~ ^ ^ . (14 

E aALd 

The time to lose the remaining energy can again be estimated by using the Larmor 
formula, with a determined from the Coulomb potential. It is easy to check that this 
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Fig. 4. Diagrams contributing to the decay B' —>■ B'-f-f. The thick 
hnes are heavy partons and the thin hnes are hght quarks. 



is dominated by the energy loss for binding energies of order E^, and we get the same 



estimate Eq. (14) for this decay time. The subsequent annihilation is much more 
rapid, as can be seen for example from the formula for the annihilation rate for an 
L = state: 

rannihilation ~ ^^""^^^0)1^ (15) 

where ip{0) is the radial wavefunction of the bound state evaluated at the origin. 
Since 'tp{0) ~ (aqcD'^)^^^ for the ground state, we have rannihilation ~ '^'^ctqcD^- We 



conclude that the decay rate for the bound state is given by Eq. (14). This decay 
occurs before nucleosynthesis (r < 1 sec) for m < 10^^ GeV, where we have used 
Ahad ~ GeV. For larger masses, the late decays to photons will affect nucleosythesis, 
and the model is ruled out |B]. 

We now discuss the decays of the bound state in the case where the heavy partons 
are electrically neutral. We assume that the brown muck is also electrically neutral. 
The energy can then be carried away by either photons or pions. The photon rate is 
a loop effect, and is suppressed by the small electromagnetic coupling, but the pion 
rate is potentially kinematically suppressed by the pion mass, which can be larger 
than the energy differences between states with AL ~ 1 for large m. We will confine 
our attention to the two photon decays, since we will see that they are sufficiently 
rapid for the most interesting range of m. 

Two photon decays occur via diagrams such as the ones in Fig. 4. Since we 
are interested in energy and momentum transfers small compared to Ahad, we can 
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parameterize the effects of diagrams such as the one in Fig. 4a by effective operators 
in the heavy parton effective theory: 

Here ip is a. heavy parton field of dimension | and is the 4-velocity of the heavy 
parton. We require two photon fields by charge conjugation invariance, and electro- 
magnetic gauge invariance then forces this to be proportional to the field strength 
tensor. Since we are interested in orbital transitions of the heavy parton, we want a 
term with at least one spatial derivative acting on the heavy parton field. The factor 
of 1/m arises because the heavy parton effective theory can depend on the 4-velocity 
and the residual momentum only in the combination = mv^ + [9j . (We ignore 
the spin indices, which are not relevant for estimating the rate for orbital transitions.) 
Note that derivatives acting on the photon fields give powers of ~ /S.E <^ Ahad, 
and so the expansion in terms of derivatives acting on the photon fields should be 
valid. Derivatives acting on the heavy parton fields give powers of l/r^, where vb 
is the size of the bound state. Since ^ 1/m, these should also be suppressed. 
There are also 4-parton processes such as the one shown in Fig. 4b that give rise to 
4-parton operators that are nonlocal on the scale -Rhad- These cannot be treated as 
local in the effective field theory we are considering, and it is more subtle to do the 



power counting for these operators. The operator Eq. ( 16 ) can be treated using stan- 
dard quantum mechanics techniques, and we are confident that it correctly counts 
the powers of m, Ahad, and in the process. We will therefore use it to get at least 
an upper bound on the lifetime of the bound state. 

Integrating over phase space, we obtain 



'■had'^B 

where is the size of the bound state. The energy of a bound state in the linear 
regime is 



E 



' ^2^2 \ 1/3 



m 



where E goes from -Emm ~ ct-Rc (where the potential becomes Coulombic, see Eq. ( 12 )) 
to -Emax ~ Ahad (where the system becomes unbound). In terms of this energy 
variable, the energy difference between AL = 1 states is 

r2 \ 1/2 

niE 



AEr^ ^\ . (19) 
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Using ~ AE, we can write the rate of energy loss as 



The hfetime is therefore 



E^FAE^ " (20) 



■dE f^B Anm^E^ 
T ^ 1 ~ / dE 



E Je. 



had 



where we have used 



EB-[^j -{AIM'/'. (22) 

Note that the hfctimc is dominated by the bound states closest to threshold, where 
AE <^ Ahad and the expansion above is valid. This decay occurs before nucleosythesis 
for m ^ 2.5 TeV, where we use Tg ~ 200 MeV and Ahad GeV to obtain the bound. 
(Note that since the rate is proportional to m^, the value of m is actually quite well 
determined.) It is possible that the decay rate to pions is more rapid, but we will not 
attempt to estimate it here. 

Let us comment on the theoretical uncertainties in our analysis above. There are 
many numerical factors that wc have estimated to be of order 1, and it is certainly 
possible that some of our estimates are numerically inaccurate. However, since we are 
interested in setting cosmological limits, the most conservative assumptions are those 
that weaken the limits. In the analysis above, wc have attempted to make "middle 
of the road" estimates for all quantities. To strengthen the cosmological limits on 
heavy stable partons, one would have to demonstrate that the estimates above are 
incorrect, taking into account the large uncertainty in hadronic quantities. Taking 
this into account, we believe that the bounds we obtain are robust. 

We now discuss the relic abundance of the heavy partons. At a temperature of 
order T ~ m/30 perturbative annihilation of heavy partons due to perturbative QCD 

gives a relic abundance Y = np/s ~ 10~^^ for temperatures T < m/30. Below 
the QCD phase transition, the second stage of annihilation described above further 
reduces the relic abundance and determines the final relic abundance. For T = Tb < 
Tc the thermally averaged rate for annihilations (more precisely, formation of bound 
states that later decay) is 

(fji ^ 1/2 
^) ' ^^^^ 
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where we expect R ~ -Rhad- This reduces the rehc abundance until the annihilation 
rate drops below the Hubble expansion rate: F = np{a\v\) < ~ gV^T"^ /M-p. 
Saturating this inequality gives a relic abundance of unbound partons of 

where s = 27r^(7*T^/45 is the entropy density, and we use ~ 10 just below the QCD 
phase transition. We are neglecting entropy production at the QCD phase transition. 

This is a very interesting relic abundance for the phenomenologically relevant 
mass range m < TeV. The cosmological bounds depend on the lifetime and decay 
modes of the heavy partons. If the lifetime is in the range of 10^ sec to 10^ sec, 
there are bounds arising from the fact that photon or hadronic decay products can 
affect nucleosynthesis [S]. For lifetimes in the range 10^ sec to 10^^ sec, there are 
bounds coming from the photodissociation of light elements [10] . For lifetimes in the 
range 10^ sec to 10^^ sec, there are bounds from distortions of the cosmic microwave 
background [llj. These bounds apply because we expect the decaying parton to have 
a significant brancing ratio into both hadrons and photons. Remarkably, all of these 
bounds become ineffective for mYp < 10~^^ GeV, near the limit of our estimated 
relic abundance for the phenomenologically relevant range of masses m < TeV. The 
conservative conclusion is therefore that such heavy partons are not excluded. For 
lifetimes in the range 10^^ sec to 10^^ sec, there are bounds from observations of 
the diffuse photon background These require that the lifetime of the partons is 
shorter than about 10^^ sec if there is a significant branching fraction into photons or 
hadrons. Stable partons are ruled out by searches for heavy hydrogen [13]. Thus, the 
reduction in the relic density by several orders of magnitude over the perturbative 
prediction lengthens the allowed lifetime for TeV scale heavy partons by roughly 12 
orders of magnitude! 
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